Toward a unified description of hadro- and photoproduction: S*— wave vr— and 

r]— photoproduction amplitudes 
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The Chew-Mandelstam parameterization, which has been used extensively in the two-body 
hadronic sector, is generalized in this exploratory study to the electromagnetic sector by simultane- 
ous fits to the TT- and 77-photoproduction S-wave multipole amplitudes for center-of-mass energies 
from the pion threshold through 1.61 GeV. We review the Chew-Mandelstam parameterization in 
detail to clarify the theoretical content of the SAID hadronic amplitude analysis and to place the pro- 
posed, generalized SAID electromagnetic amplitudes in the context of earlier employed parameterized 
forms. The parameterization is unitary at the two-body level, employing four hadronic channels and 
the 7A'^ electromagnetic channel. We compare the resulting fit to the maid parameterization and 
find qualitative agreement though, numerically, the solution is somewhat different. Applications of 
the extended parameterization to global fits of the photoproduction data and to global fits of the 
combined hadronic and photoproduction data are discussed. 

PACS numbers: 13.60.-r, 11.55.Bq, ll.80.Et, ll.80.Gw, 13.60.Le 



I. INTRODUCTION 

Most of our knowledge of the excited baryons has come 
from fits to hadronic scattering data, in particular pion- 
nucleon scattering, for which there exists an accurate and 
nearly complete database extending through and above 
the resonance region. Sufficient polarization observables 
exist to constitute complete measurements over a sig- 
nificant kinematic interval. The range of hadroproduc- 
tion data, including TriV — > ttN, ttN — )■ rjN, ttN 
ljN , and other inelastic processes including, for example, 
strangeness production, have also been used to constrain 
theoretical models and phenomenological parameteriza- 
tions of the scattering and reaction amplitudes. 

Currently, however, a renaissance is underway in me- 
son production and resonance physics with reaction data 
issuing from a number of precision electromagnetic fa- 
cilities. Collaborative theoretical and phenomenological 
efforts have started to analyze these data in ways consis- 
tent with some subset of constraints imposed by quantum 
field theory upon the reaction amplitudes. The quality 
and quantity of data in electromagnetic induced reactions 
is becoming sufficient to rival and possibly surpass the 
hadroproduction data. Since the electromagnetic reac- 
tions proceed mainly through the hadronic channels, the 
new data offers the possibility of "back-constraining" the 
hadronic amplitudes, conventionally determined only in 
fits to the hadroproduction data. 

It is in this context that we have completed an ex- 
ploratory study of the S'-wave tt- and ry-photoproduction 
multipoles in the "Chew-Mandelstam" approach, related 
to the N/D representation, to the electromagnetic reac- 
tion amplitude. The novel concept, which provokes and 
permits this exploratory study, is the generalization of 
the Chew-Mandelstam approach to the electromagnetic 



sector. We have developed a new form for the amplitude 
that incorporates multichannel hadronic rescattering ef- 
fects in a complete manner consistent with unitarity. The 
near-term objective is to develop a framework in which 
to analyze the hadro- and elect reproduction reactions si- 
multaneously in a global framework. 

Recent experimental observations of the photoproduc- 
tion of the rj meson from the proton have yielded mea- 
surements of the unpolarized differential cross section ^ 
3 and photon beam asymmetry[4j _5j of high precision. 
Forthcoming measurements from the CLAS Collabora- 
tion at Jefferson Lab [5] and Mainz [7] will rival, if not 
surpass, the precision of the existing measurements. 

Several interesting features of r] meson physics moti- 
vate these measurements and their theoretical interpre- 
tation in various fields of nuclear physics, astrophysics, 
and particle physics. The possibility that the 77-nucleon 
interaction may be attractive milj suggests the existence 
of bound states of the 77 meson with nuclei. Certain res- 
onances, the 5*11(1535) A'^* resonance in particular, are 
significantly coupled to the rjN channel, and the pho- 
toproduction of this final state provides an independent 
method to probe the isospin T = ^ resonance spectrum 
and its couplings [10] . 

The strong interactions of the tt and ry mesons require 
multichannel descriptions which respect unitarity in the 
relevant channel space in order to obtain a realistic de- 
scription of the data. The Chew-Mandelstam iC-matrix 
approach [TTI - [T3] is an effective parameterization of the 
observed reaction data since the elements of the Chew- 
Mandelstam K matrix may be assumed to be real if the 
couplings to neglected, open channels are small. 

Several relatively recent K matrix analyses of the cou- 
pled ttA^, rjN , and channels |14f[r7] have been suc- 
cessful in obtaining reasonable parameterizations of the 
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two-body partial wave amplitudes [TS]. The purpose of 
the present work is to investigate the extent to which a 

1/2 

description of the tt photoproduction Eq_^_ {Sh) ampli- 
tude and the modulus of the rj photoproduction ampli- 
tude yields an 77 photoproduction multipole with a res- 
onant phase. Various calculations [El fT9H2T] . indicate 
that the modulus of the rj photoproduction amplitude 
near threshold is fairly model-independent, being repro- 
duced in a range of calculational models or schemes. In 
the present work, we take as input hadronic T matrix ele- 
ments, determined in realistic (^-squared per datum ~ 1) 
fits to data [22], as discussed in the following sections. 

In Sec. [n] we review in some detail the Chew- 
Mandelstam form|13j of the parameterization. The pur- 
pose of this review is to establish the theoretical consid- 
erations that motivate the amplitude parameterizations 
used in the SAID program, to place these amplitudes in 
the context of other hadronic amplitude parameterization 
schemes and to lay the groundwork for future improve- 
ments. Section |III| gives the results for the fits to the 
isospin T — ^ tt -photoproduction amplitude, -EJ^, and 
the modulus of the rj photoproduction amplitude, i?o+- 
The conclusions are given in Sec. |IV| We find, in this ex- 
ploratory study, an 77 photoproduction multipole having 
a resonant shape, qualitatively similar to a Breit-Wigner 
form, and similar to other calculations [19 1 [211 123j . There 
is, however, significant deviation from the simple Breit- 
Wigner form. 



II. CHEW-MANDELSTAM 
PARAMETERIZATION 

Previous work in the determination of the i] photopro- 
duction amplitudes |14l[l6] has shown that an approach 
which includes the coupling of the electromagnetic chan- 
nel to the TriV and t]N channels in the region of energies 
near the center-of-mass energy, W = 1535 MeV gives a 
reasonably good description of the data and a plausible 
form for the amplitudes. However, as our ultimate objec- 
tive is the simultaneous parameterization of hadro- and 
photoproduction scattering and reaction observables, we 
will go beyond the two-channel treatment for this study 
of the Eq, multipole amplitude. 



A. Unitarity constraint 

The form of the Chew-Mandelstam parameterization, 
which we employ in this study follows as a consequence 
of the analytic structure imposed by the unitaritv [21H25] 
of the S matrix in the physical region, W > uii + mj, 
where W is the center-of-mass energy and and mt are 
the masses of the incident and target particles. Confining 
our attention to two-particle initial and final states, the 



S matrix is defined as 

Sc^piE) = (k„a|5|k^/3) (1) 

+ 2t7:6{E^-Ef3)(k^a\T\kp(3) (2) 

where k^^^ are the final and initial relative momenta, 
respectively, E — Ea — Ep — W is the center-of-mass 
energy, and the labels a and /? denote the particle species, 
spins, and internal quantum numbers, such as isospin. 
The initial and final energies, Ep and E^, respectively 
are related to the on-shell relative momenta for channel 
, Q as 

W = + S„,2 (3) 
Ec.,^ = ,J~kl+ml,. (4) 

The on-shell relative momentum may be expressed in 
terms of the center-of-mass energy, W, as 

ka = - ma+y^W - 

X v^W^ + m^+^iy + m„_, (5) 

with ma± — nia.! ± ma,2- 

The scattering operator, S, is unitary 

S^S ^SS^ ^1 (6) 

and if we restrict our analysis to energies where just two- 
particle channels contribute, we obtain 

/d'fca(k„a|5^|k „a){'kaa\S\ki3/3) 

^S'^^\ko,~kp)S^p. (7) 

Substitution of Eq. ([2| into the relation above yields the 
unitarity constraint on T 

Ta.p - Tip = 2TnJ2 [ d\„Tij{E^ - E,)T„p. (8) 
Effecting the integration on fco- = |ko-| gives 

-T1p = 2^Y. [ d^AHW - m,+)p,T,0 (9) 

where 

P'y(ka) = • (10) 

The presence of the Heaviside step function, 9{W — mcT+) 
is a consequence of the fact that, over the range of inte- 
gration ka- > 0, the argument of the 6 function has a 
solution, Ea — E^ = only when W > ma-+. Equation 
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^ implies discontinuities in the derivative of the imagi- 
nary part at each channel threshold W = ma+'- 



Y- i^ais - T*p] - Im TafS 



(11) 



E 



Jdn,T*JiW - m,+ )p„T,p (12) 



where we have assumed that, due to the time-reversal 
invariance of the strong interaction, T^p = Tp^- The 
violation of the Cauchy-Riemann equations at threshold 
indicates the presence of a branch point. We distinguish 
between the dynamical singularities at each threshold 
opening mg-^ and kinematical singularities, due to the 
presence of kinematical factors such as k^. The kine- 
matical singularities are removed from the unitarity con- 
straint by considering T'^^ = ^^T^^^^. 

We may transform to the partial wave representation 
and write 

r^p ~ = 2z ^ T':J{W - m,+ )T',p (13) 



where the T^^ now represent the partial wave amplitudes. 
Casting this relation as a matrix equation 



. [T' ~ T'*] = T'*e{W - M+)T', 



(14) 



where M^ ^a- = 'm-rr+^atT, s^nd multiplying from the left 
by [T'*]~^ and from the right by T'^^ gives 



Im T' 



-e{w - A'U 



(15) 



a diagonal matrix. Since this equation isolates the imag- 
inary part of the inverse-T matrix, we write 



T'-^ = Rc T'^i +ilm T'"\ 



(16) 
(17) 



where we've defined Re T' ^ 



K' 



-1 and ^ 
^ PaKafj^J Pp. Multiplying from one side by T' and the 
other by K' gives the Heitler integral equation |291 15(7| 



T' ^K' + K'i6{W - M+)T'. 



(18) 



This is the starting point for the Chew-Mandelstam pa- 
rameterization of the reaction amplitude. 

We emphasize that, in the physical region, the uni- 
tarity relation is satisfied by the imaginary part of 
T'^^ . Therefore the Heitler K matrix is analytic, except 
for possible isolated poles [ST, throughout the physical 
region [Ml I32j . This is apparent if we consider a dynam- 
ical equation of, for example, the Lippmann-Schwinger 
form: 



T = V + VGqT, 
h, — Uq 



(19) 
(20) 



Here V is the interaction part of the full Hamiltonian, 
E = VF, Hq is the free-particle Hamiltonian, and V de- 
notes the Cauchy principal value prescription. Substitu- 
tion of Eq.([20| into Eq.([l9| gives T = K + iK5{E-Ha)T 
where 



K ^V + VV 



1 



E-Ho 



-K. 



(21) 



The Cauchy principal value prescription in this equation 
yields a kernel which is completely continuous in the 
physical region. The spectrum of the kernel therefore 
possesses no eigenvalues in the continuum and K is ana- 
lytic (other than possible poles) there [331 • The K matrix 
may possess singularities in other regions of the complex 
energy plane. In fact, the interaction V possesses singu- 
larities in regions outside the physical region In 
particular, there is a branch point at some value W < Q. 
We intend to neglect singularities in the region Re < 
for the purposes of the present study and avoid a detailed 
discussion of them here. A description is available in the 
literature j34H37) . Inclusion of singularities in the region 
Re < will be explored in subsequent investigations. 

The partial wave amplitude is therefore known to have 
the following singularities. There are branch points in 
the physical region at the channel-opening thresholds as 
in Eq.Q, branch points in the region W < and pos- 
sible poles consistent with causality [32l [38]. An efficient 
parameterization following Ref. [37j . which encodes these 
singularities, involves the factorization of the partial wave 
amplitude. This is referred to as the "N/D" approach. 
We will use the N/D language to clarify the nature of 
the singularities of the T matrix which are included and 
those neglected in our Chew-Mandelstam approach. 



B. Relation to N/D approach 

The N/D approach has been used to analyze a variety 
of reactions[371 [321 HO]- As our long term objective is 
the generalization of the existing method used to param- 
eterize the hadronic and electromagnetic amplitudes, we 
collect here some of the relevant equations of the N/D 
approach. The T matrix is written in the factorized form 



T{W) = D-'^(W)N{W) 



(22) 



where N and D are Nch x N^h arrays [28 , where N^h is 
the number of included two-body channels. This relation 
has been shown to be consistent with the requirement of 
time-reversal invariance in Ref. [41]. The relations 

Im D{W) = N{W)lmT-\W) W > m, + rat (23) 
ImiV(iy) = W >m,^mt (24) 

\mN(y^) = D{y/)\mT{W) iy<0 (25) 



\mD(y^) = 



ly < (26) 



give the essential content of the N/D approach. They 
state that the function D has branch points only in the 
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physical, W > rrii + rrit region and that A'' has only un- 
physical, < branch points. These relations deter- 
mine the following dispersion relation (or Hilbert trans- 
form) representation for D 



rir, 



D{w) = - 1 - 

1=1 



TT 



oo 



N{W')p{W') 



(27) 



with Up subtractions. Here, D{W; Wi) is a polynomial 
of order rip, W G C, and Wt is the lowest production 
threshold. Here, we show the polynomial ambiguity of 
the Hilbert transform explicitly to allow for the possibil- 
ity that the parameterization includes several subtraction 
points. 

Using the relation T = ND^^ in the physical region, 
the numerator factor, N can be shown to satisfy the in- 
tegral equation 



N{W) 



\ 1 



1=1 



N{W')p{W') 



(28) 



where f denotes the Cauchy principal value integral, and 
the Heitler K matrix, K is defined by Eq.( 18 1, with K ~ 

-^1 T^l -1 

C. Chew-Mandelstam parameterization 

The preceding discussion of unitarity and the N/ D ap- 
proach provides the context for our present parameteri- 
zation. The Chew-Mandelstam parameterization devel- 
oped here is similar to those of Refs. |lllfT^ and [T3]. We 
consider Eq.(17) and rewrite it, confining our attention 



to the S wave multipole as 



T-^ = R-^ -ip 



{K- 



f ReC) 



(Re C + ip) 



(29) 



(30) 



where p = pd{W - M+) and ImC = p = e{W - M+)p. 
The transition matrix is given in terms of the "Chew- 
Mandelstam" (CM) K matrix, K by 



T = K + KCT. 



(31) 



Equation (31 1 fixes our Chew-Mandelstam parameter- 
ization. In the language of the preceding section, we 
have neglected the W < branch points of N and 
made the approximation N{W) = K{W), an entire func- 
tion. The "Chew-Mandelstam" function, Cq is deter- 



Chew-Mandelstam function is given by a Cauchy integral 
over the discontinuity of Cq in the physical region 

dW Pa{W') 



^ W -W 

°° dW Pa{W') 



Wt 



TT W'-W, 



(32) 



where we have made one subtraction, < Wg < Wt- 
Defining Za 



w-w ^® '"^^ rewrite Eq.(32) as 
''^ dx p{x) 



Jo 7^ 



Zc.(M^)' 



(33) 



The relationship between the Heitler K matrix and the 
CM K matrix, K is given by 



K = K + K[ReC]K. 



(34) 



This demonstrates a possible advantage of using the CM 
K matrix. If we consider a polynomial parameterization 
of a given CM K matrix element 



K 



E-n 

n=0 



(35) 



where n^^ are channel dependent integers controlling the 
order of the polynomial (polynomials typically less than 
fifth order are used) and z^p is a possibly channel de- 
pendent linear function of the center-of-mass energy, W 
then we see, by solving Eq.(34) for K 



K 



' -K 



l^K[ReC] 



(36) 



mined solely by the unitarity constraint, Eq.(15l since 
Eq.(30l is equivalent to taking D = 1 — KC. Then the polynomial approximation for K. 



that poles may appear in the K matrix. Attempts 
to relate the K matrix poles to resonances have been 
made P^Hiij . Here, we simply point out that, though K 
matrix poles are not simply related to T matrix poles[3T], 
Eq.(p4|) shows that one need not explicitly include pole 
terms in K in order to have poles in K. Parameteriz- 
ing K{W) — N[W) as a polynomial, as noted, neglects 
singularities in the unphysical region, W < 0|33]- The 
branch points there and discontinuities across their as- 
sociated branch cuts are determined by the production 
mechanisms 46J relevant for the reaction considered. 

There are at least two reasons why polynomials may 
provide a reasonable starting point for a realistic param- 
eterization of multichannel scattering and reaction am- 
plitudes. The unitarity branch points, given their phys- 
ical nature, largely determine the gross structure of the 
amplitudes in the physical region. This leads in an ob- 
vious way to the supposition that more distant singu- 
larities in the complex W plane associated, in partic- 
ular, with the branch points in the unphysical region 
may be less important. Experience has also confirmed 
this to be true. The existing said parameterizations of 
ttN elastic scattering[55], the TriV — 7> jyA^ITU] reaction, 
TT-photoproduction j47l l48] and electroproduction[49] and 
other reactions all reveal that a realistic description with 
per datum in the range of 1 — 3 is possible with the 
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SP06 


FA02 


KA84 








2.0 


6.1 


2.1 


8.8 


5.0 


24.9 


1 




'p - 


-¥ TT^p 


1.9 


6.2 


2.0 


6.6 


9.1 


51.9 






'p - 




2.0 


4.0 


1.9 


5.9 


4.4 


8.8 




TV' 


'p - 


rjn 


2.5 


9.6 


2.5 


10.5 









TABLE I: Normalized (left of each column pair) and un- 
normalized (right of each column pair) x'^'P^r-datum for 
the SP06[22] and FA02[50] solutions of SAID, KA84[5T]. 
EBAC[52], and Giessen[55]. The energy ranges of the four 
groups are from threshold to 2.5, 2.9, 1.91, and 2.0 GeV, 
respectively. [54] 



III. RESULTS 

The Chew-Mandelstam parameterization for the T ma- 
trix, described in the preceding section, has been applied 
recently [22] to a coupled-channel fit for the ttN elastic 
scattering and TriV rjN reaction. It gives a realistic 
description of the data with P^r datum better than 
any other parameterization or model, to our knowledge. 
The psr datum is shown in Table ^ against other pa- 
rameterizations and model calculations for which we pos- 
sess sufficient amplitude information to perform such an 
analysis [54|. The current said parameterization used in 
this fit is given as 



al3 



(37) 



where a, /3 and a are channel indices for the considered 
channels, ttNjTtA, pN and i]N. This parameterization 
has been discussed in Refs.|131 [Til [55] . Given the suc- 
cess of this approach in the hadronic two-body sector, 
the application to the study of meson photoproduction 
is warranted. 

The central result of the current exploratory study is 
to show that this form can be extended to include the 
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FIG. 1: The SAID 5*11 multipole for the ttN — > tjN reaction 
as a function of energy, W' lQ., . The solid (dashed) line is the 
real (imaginary) part of the amplitude. 
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FIG. 2: Comparison of the SAID real (solid curve) and imagi- 
nary (short-dashed curve) parts of the Eq^ multipole ampli- 
tude with that of the MAID [55] real (long-dashed curve) and 
imaginary (dot-dashed curve) parts. The amplitudes are plot- 
ted along with the real (circles) and imaginary (squares) of 
the SAID single-energy solutions |22). 
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FIG. 3: The predicted values for the real (solid curve) and 
imaginary (dashed curve) for _EQ^versus the energy, W. The 
modulus I-Eq+I (dotted curve), the real (dot-dashed curve), 
and the imaginary (double dot-dashed curve) parts of the tt- 
photoproduction, Eq^ were fit to pseudodata generated from 
the SAID solution [47] with the parameterized form Eq.|38l 
using 8 parameters (see text). 



electromagnetic channel. 



(38) 



where 7 denotes the electromagnetic channel, jN. Note 
that Eqs.(37) and (38) share the common factor, [1 — 
KC]~^ which encodes, at least qualitatively speaking, 
the hadronic channel coupling (or rescattering) effects. 

The form Eq.(38l for photoproduction should be 
contrasted with that currently employed in the tt- 
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FIG. 4: The r;— photoproduction Sn multipole amplitude, 
-Eq^ versus the energy, W fit using the previously employed, 
non-unitary form of Eq.(39|. The behavior near W ~ 1535 
MeV is not resonant as can be clearly seen in Fig.Q. 



photoproduction studies of Refs. |151 HTl H5] 

T^^ = A{W){1 + tT^AW)) + iB{W)T^jW) 



(39) 



are 



where the "structure functions" A(M^) and B{W) 
parameterized as polynomials in the energy, W, T^^^ — 
T-„N,-yN and Ttttt = Tttn.tiN, and the factor yl(VF) contains 
a contribution from tree-level Born diagrams. This sat- 
isfies Watson's theorem 57 (as does Eq.(38)), and is de- 
rived via the considerations discussed in Ref. [58 . While 
resulting in a realistic description of the data and being 
comparable, at least qualitatively, with other parameter- 
izations such as MAID [23] for tt -photoproduction, it does 
not satisfy the full multichannel unitarity constraint im- 
posed by Eqj9]). This deficiency led us to consider the 
form in Eq. ( |38[ ), which manifestly satisfies the multichan- 
nel unitarity constraint, Eq.([9|. 

The need to include the multichannel unitarity effects 
of Eq. ( 38 1 have also become apparent in difficulties faced 



in attempts to parameterize the 77-photoproduction re- 
action using forms 59J similar to Eq.(39). Forms of this 



type, used in fits to the ry-photoproduction data alone, 
yielded an S'-wave multipole without a clearly resonant 
shape, even while yielding fits to the observed data with 
realistic per datum on the order of 2 to 4. An example 
of such a fit employing Eq.(39) is shown in Fig.Q. Near 
values of the center-of-mass energy W ~ 1535 the am- 
plitude in Fig.Q is decidedly not resonant. This is also 
clear in the Argand plot of Fig. ([5]). Here we have shown 
the comparison of the fit forms used in Ref . ,59J (with en- 
ergies marked by triangles) This difficulty was an early 
motivation for the present study. Expectation of reso- 
nant behavior for 77-photoproduction, — >■ rjN in the 
S'-wave can be argued straightforwardly. For example, 
since the electromagnetic coupling to the ttN channel 
is large, the 7^'^ — >■ 'qN reaction may proceed via the 
ttN — >■ rjN amplitude of Fig.Q, or through direct reso- 
nance production. We therefore anticipate the hadronic 



subprocess will drive a significant resonant effect in the 
(isoscalar) electromagnetic transition. 

Several other studies have determined that a res- 
onant structure near W ~ 1535 MeV is consistent 
with the reaction data. These works include those in 
Refs.[Tni m HD We should note that all of these 
works have assumed the Sn wave to be resonant, usu- 
ally by including a Brcit-Wigner or similar term explicitly 
into their formalism. We do not make this assumption 



in using Eq.(38 I 



In light of the study of Ref. 59J and the necessity of in- 
cluding the full multichannel unitarity for the purposes 
of obtaining a global description of the hadro- and pho- 
toproduction data, we have carried out an exploratory 
study to determine the efficacy of doing such a fit within 
the Chew-Mandelstam parameterization, Eq. p8| . In the 
present study we perform a coupled-channel fit of the 
modulus |£^Q_^(VF)| and (the real and imaginary parts of) 
the existing said and maid 7r-photoproduction ampli- 
tudes, E^^ in the S'ii(1535) resonance region, compared 
in Fig.(l2j. The fit was carried out by taking the factor 
[1 — KC]'^^ in Eq.(38) as determined in the the hadronic 
study of Ref.m] and adjusting the parameters of if 0-7 
(discussed in detail in the subsections below) . The phase 
of the EQ_^_ multipole in this study gives a resonant wave 
and encourages us to continue with this approach, as dis- 
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FIG. 5: Argand-plot comparison of the 77— photoproduction 
Sii multipole amplitudes, HeEg^ versus ImEg^ plotted in 
the range 1490 MeV < W < 1610 MeV of center-of-mass 
energy, W with two fit forms. The curve with energies 
marked by triangles is another representation of the result 
for _Eq_,_ shown in Fig.Q, determined using the parameteriza- 
tion of Eq.( |39[ ). The curve with energies marked by circles is 
another representation of the result for Eq_^_, shown in Fig.([3]l, 
determined using the parameterization of Eq. ( 38 1 . The curves 
span the same interval in energy but with different spacings. 
The first curve (triangles) is clearly non-resonant in the region 
shown while the second curve (cirlces) clearly shows resonant 
behavior; the apex on the Argand diagram of the second curve 
occurs at precisely W = 1535 MeV. 
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FIG. 6: The predicted values for the real (solid curve) and 
imaginary (dashed curve) for E'p.^versus the energy, W. The 
modulus |i?o+l (dotted curve), the real (dot-dashed curve), 
and the imaginary (double dot-dashed curve) parts of the tt- 
photoproduction, _Eo+ were fit to pseudodata generated from 
the MAID solution |60| with the parameterized form Eq.(|38[) 
using 7 parameters (see text). 



FIG. 8: The real (solid curves) and imaginary (dashed curve) 
for i?g_|_from the 7-parameter fit in Fig.(|6| compared with the 
t;-maid solution [23], marked by squares. 



cussed in the final section. 

The decision to fit the modulus |i?o+l based on em- 
pirical considerations. The MAiD[23j parameterization by 
the Dubna-Mainz- Taipei Collaboration (DMT) and the 
model calculations in Refs.UHl EHl HH HI] agree at the 
few-percent level on the modulus of the low-energy r/~ 
photoproduction amplitude, |£'Q^(iy)|. This is antici- 
pated on the grounds that, in the 511(1535) resonance 
region, the differential cross section is largely angle in- 




T490 1510 1530 1550 1570 1590 1610 
W [MeV] 

FIG. 7: The predicted values for the real (solid curve) and 
imaginary (dashed curve) for i?Q^versus the energy, W. The 
modulus \Eg_^_\ (dotted curve), the real (dot-dashed curve), 
and the imaginary (double dot-dashed curve) parts of the tt- 
photoproduction, Eq^ were fit to pseudodata generated from 
the MAID solution [60] with the parameterized form Eq.(38l 
using 14 parameters (see text). 



dependent and therefore dominated by the S'-wave pro- 
duction. It also indicates that the production is largely 
resonant, but we do not make this common assumption. 

While the modulus \Eq_^_\ appears to be known at 
the level of a few percent, the 7r-photoproduction 5*11 
amplitude is, surprisingly, not very well determined 
through different parameterizations. Figure ([2| shows 
the SAID [22] and MAID [23] results for Eq^. Given this 
discrepancy, we have also carried out the fit described 
above with the modulus \Eq_^_\ and the maid parameter- 
ization. 




1530 1550 1570 
W [MeV] 



1610 



FIG. 9: The real (solid curves) and imaginary (dashed curve) 
for iSg^from the 14-parameter fit in Fig.([7| compared with 
the 77-MAlD solution[23]. marked by squares. 
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A. Fit with SAID 

Figure ([s]) shows the resuh of fitting the modulus 
\Eq^\ and the real and imaginary parts of the SAID Eq^ 
multipole 07] using an eight-parameter fit. The Chew- 
Mandelstam K matrix was assumed to have the form 



(40) 



taking n^^ = 1, in Eq.(35l, for a and /3 taking values 



in the set of four channels, ttN, ttA, pN, and ijN. The 
energy variable Zap is 



Zat3 = W -Wt,a, 



(41) 



where the threshold masses, Wt^a are + rriM, Sm^r + 
rriN, 2m^ +mN, and + toat for a = ttN, ttA, pN, and 
TjN, respectively, and Wt^a is taken to be the lower of the 
thresholds for channels a and /3. The eight parameters 
were varied in the fit to a total of 113 pseudodata points 
including the modulus I-Bq+I over the energy range 1490 
MeV <W < 1610 MeV and the amplitude E^_^ over the 
energy range 1120 MeV < W <1610. The per datum 
over for the fits to the pseudodata, generated with the 
SAID interactive code facility were less than one in all 
of the fits made in this work including those in the region 
1120 MeV <W < 1490 MeV which are not displayed in 
order to keep the figures manageable and focus attention 
on the 511(1535) resonance region. The pseudodata were 
assigned 5% errors in the fit. 



B. Fit with MAID Eq^ 



The graphs in Figs. ([6]) and ([7| used seven and four- 
teen parameters, respectively, to fit the \Eq^\ and 
i?o_,_ amplitudes from MAID [23]. The seven-parameter fit 
in Fig.(|6| is the minimal set of parameters needed to 
obtain a per datum < 1. The parameters used in 



Eq.(35l for this fit were c,rj,n, n = 0,1,2, Cp^,„, where 
n = 0, 1 and Cjjj^ and Crj-f^i. The quality is degraded at 
the higher energy end of the fit region for the imaginary 
part of Eq_^. Nearly perfect agreement is obtained if we 
use a fourteen parameter form for K^r-y ■ The parameters 
used in Eq.(35) for this fit were c^^^n, ca7,ti, Cp-y^„, where 



0; 1, 2, 3 and c^-y and c™^i 



Note, from Figs.([8|) and (|9|, that the fit giving the 
better representation of the maid amplitude is sim- 
ilarly closer to the maid i?Q_|_ result. This is somewhat 
surprising perhaps, since although the maid pion- and 
eta-photoproduction use the same pole positions in both 
amplitudes, these parameterizations arc not constrained 
by unitarity. 



IV. CONCLUSION AND ONGOING WORK 

We reviewed the implication of unitarity on the an- 
alytic structure of the single meson production scatter- 



ing and reaction amplitudes. The non-analyticities in 
the regions W > and < 0, the right- and left- 
hand cuts, respectively were demonstrated to be prop- 
erly accounted for by the N/D approach. We related 
the Chew-Mandelstam iiT-matrix parameterization to the 
N/ D approach, showing that the parameterization of the 
K matrix neglects the effects of the distant left-hand 
cut. The purpose of this review is to place the long-used 
SAID amplitudes in the context of other hadronic ampli- 
tude parameterization schemes and to lay the ground- 
work for future improvements to the existing parameter- 
ization forms. 



Using the Chew-Mandelstam K matrix K, we per- 
formed a simultaneous coupled-channel fit of the 77- 
photoproduction multipole modulus, I-Eq+I ^^'^ 
TT-photoproduction amplitude, EQ_^_. The parameteriza- 
tion was restricted only to the CM K matrix elements 
Kcrj in Eq.(38), while the [1 — KC]^^ factors were taken 
from the existing SAID fits to the hadronic data. The an- 
ticipated resonant structure for the phase of the Eq^ mul- 
tipole was demonstrated in fits to both said and maid 
amplitudes. 

The results of the exploratory study indicate that this 
is a reasonable approach toward the objective of deter- 
mining a complete set of scattering and reaction am- 
plitudes for ttN — > TriV, ttN — )• r]N, — )• TriV, and 
7iV rjN processes in a multichannel unitary formal- 
ism. The first stage in this procedure, demonstrating 
that coupled-channel simultaneous fits of the tt- and t]- 
photoproduction reactions for a single partial wave (6*11) 
is possible, has been completed. The next phase, consists 
of a fit to the tt- photoproduction reaction observables. 
Following this, a simultaneous fit to the reaction observ- 
ables for the TT- and ry-photoproduction reactions will be 
performed. As a practical matter, these two phases will 
be completed using the [1 — i^C]^^ "rescattering" factors 
determined in separate fits to the hadronic scattering and 
reaction data. The final phase of the study will be a si- 
multaneous fit to both the hadronic and electromagnetic 
scattering and reaction observables and will constitute, 
at least for two-body unitarity, a global description of the 
hadro- and photoproduction amplitudes. 



Acknowledgments 



The authors thank R. Arndt without whom this work 
would not be possible. This work was supported in 
part by the U.S. Department of Energy Grant DE-FGG2- 
99ER41110. We thank the Department of Energy's Insti- 
tute for Nuclear Theory at the University of Washington 
for its hospitality and the Department of Energy for par- 
tial support during the initiation of this work. 



9 



[1] V. Crede et al. (CB-ELSA), Phys. Rev. Lett. 94, 012004 
(2005), hep-ex/0311045. 

T. Nakabayashi et al., Phys. Rev. C74, 035202 (2006). 
M. Williams et al. (CLAS), Phys. Rev. C80, 045213 
(2009), 0909.0616. 

J. Ajaka et al., Phys. Rev. Lett. 81, 1797 (1998). 
D. Eisner et al. (CBELSA), Eur. Phys. J. A33, 147 
(2007), nucl-ex/0702032. 
_M. Dugger et al. (CLAS) (2004), URL|http : //www . j lab^| 
org/ exp_prog/proposals/05/PR05-012 .ps 

C. McGeorge et al., Eur. Phys. J. A37, 129 (2008), 
0711.3443. 

W. K. Cheng, T. T. S. Kuo, and G. L. Li, Phys. Lett. 
B195, 515 (1987). 

Q. Haider and L. C. Liu, Acta Phys. Polon. Supp. 2, 121 
(2009), 0902.4248. 

S. Prakhov et al., Phys. Rev. C72, 015203 (2005). 
J. L. Basdevant and E. L. Berger, Phys. Rev. D19, 239 

(1979) . 

B. J. Edwards and G. H. Thomas, Phys. Rev. D22, 2772 

(1980) . 

R. A. Arndt, J. M. Ford, and L. D. Roper, Phys. Rev. 
D32, 1085 (1985). 

A. M. Green and S. Wycech, Phys. Rev. C55, 2167 
(1997), nucl-th/9703009. 

R. A. Arndt, A. M. Green, R. L. Workman, and 
S. Wycech, Phys. Rev. C58, 3636 (1998), nucl- 
th/9807009. 

A. M. Green and S. Wycech, Phys. Rev. C60, 035208 
(1999), nucl-th/9905011. 

V. Crede et al. (CBELSA/TAPS), Phys. Rev. C80, 
055202 (2009), 0909.1248. 

The K matrix approach of Refs. |14t416) include contribu- 
tions to the hadronic T matrix which represent the ttttN 
channel, though the three-body cut is not accounted for 
properly. 

N. Kaiser, T. Waas, and W. Weise, Nucl. Phys. A612, 
297 (1997), hep-ph/9607459. 

L. Tiator, G. Knochlein, and C. Bennhold, PiN Newslett. 
14, 70 (1998), nucl-th/9802064. 

L G. Aznauryan, Phys. Rev. C68, 065204 (2003), nucl- 
th/0306079. 

R. A. Arndt, W. J. Briscoe, I. L Strakovsky, and 
R. L. Workman, Phys. Rev. C74, 045205 (2006), nucl- 
th/0605082. 

W.-T. Chiang, S.-N. Yang, L. Tiator, and D. Drechsel, 
Nucl. Phys. A700, 429 (2002), nucl-th/0110034. 
W. Zimmerman, Nuovo Cim. 21, 249 (1961). 
R. J. Eden, Proc. Royal Soc. of London 210, 388 (1952). 
J. C. Polkinghorne, Nuovo Cim. 23, 360 (1962). 
J. B. Boyling, Nuovo Cim. 33, 1356 (1964). 
J. D. Bjorken, Phys. Rev. Lett. 4, 473 (1960). 
W. Heitler, Math. Proc. Camb. Phil. Soc. 37, 291 (1941). 
M. L. Goldberger and K. M. Watson, Collision Theory 
(John Wiley and Sons, Inc., New York, 1964). 
R. L. Workman, R. A. Arndt, and M. W. Paris, Phys. 
Rev. C79, 038201 (2009), 0808.2176. 
R. J. Eden, P. V. Landshoff, D. L Olive, and J. C. Polk- 
inghorne, The Analytic S-matrix (Cambridge University 
Press, Cambridge, 1966), pp. 231-232. 
S. Weinberg, Phys. Rev. 133, B232 (1964). 



[34] S. Mandelstam, Phys. Rev. 112, 1344 (1958). 
[35] S. C. Frautschi and J. D. Walecka, Phys. Rev. 120, 1486 
(1960). 

[36] W. R. Frazer and J. R. Fulco, Phys. Rev. 119, 1420 
(1960). 

[37] G. F. Chew and S. Mandelstam, Phys. Rev. 119, 467 

(1960) . 

[38] N. M. Queen and G. Violini, Dispersion theory in high- 
energy physics (John Wiley & Sons, Inc., New York, 
1974), pp. 13-14. 

[39] R. A. Arndt, Phys. Rev. 165, 1834 (1968). 

[40] A. O. Barut, The theory of the scattering matrix (The 
Macmillan Company, New York, 1967), pp. 211-222. 

[41] J. D. Bjorken and M. Nauenberg, Phys. Rev. 121, 1250 

(1961) . 

[42] A. D. Martin and T. D. Spearman, Elementary parti- 
cle theory (North Holland Publishing Co., Amsterdam, 
1970). 

[43] R. M. Davidson and N. C. Mukhopadhyay, Phys. Rev. 

D42, 20 (1990). 
[44] S. Ceci, A. Svarc, B. Zauner, M. Manley, and S. Capstick, 

Phys. Lett. B659, 228 (2008), hep-ph/0611094. 
[45] We note that the current parameterization neglects the 

constraints of chiral symmetry near threshold and in the 

unphysical region. Therefore the Adler zero[63] is, too, 

neglected. 

[46] O. Babelon, J. L. Basdevant, D. Caillerie, and G. Men- 
nessier, Nucl. Phys. B113, 445 (1976). 

[47] R. A. Arndt, R. L. Workman, Z. Li, and L. D. Roper, 
Phys. Rev. C42, 1853 (1990). 

[48] R. A. Arndt, R. L. Workman, Z. Li, and L. D. Roper, 
Phys. Rev. C42, 1864 (1990). 

[49] R. A. Arndt, I. I. Strakovsky, and R. L. Workman, PiN 
Newslett. 16, 150 (2002), nucl-th/OUOOOl. 

[50] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, R. L. Work- 
man, and M. M. Pavan, Phys. Rev. C69, 035213 (2004), 
nucl-th/0311089. 

[51] R. Koch, Z. Phys. C29, 597 (1985). 

[52] B. Julia-Diaz, T. S. H. Lee, A. Matsuyama, and T. Sato, 
Phys. Rev. C76, 065201 (2007), arXiv:0704.1615 [nucl- 
th]. 

[53] V. Shklyar, H. Lenske, U. Mosel, and G. Penner, Phys. 
Rev. C71, 055206 (2005), nucl-th/0412029. 

[54] The per datum analysis presented in Table (|l| in- 
cludes results for groups that have provided amplitudes 
in a manner allowing the reconstruction of observables at 
arbitary angle and energy. 

[55] D. Drechsel, S. S. Kamalov, and L. Tiator, Eur. Phys. J. 
A34, 69 (2007), 0710.0306. 

[56] R. A. Arndt, I. I. Strakovsky, R. L. Workman, and M. M. 
Pavan, Phys. Rev. C52, 2120 (1995), nucl-th/9505040. 

[57] K. M. Watson, Phys. Rev. 88, 1163 (1952). 

[58] R. L. Workman, Phys. Rev. C74, 055207 (2006), nucl- 
th/0510025. 

[59] R. A. Arndt, private communication (2009). 

[60] W.-T. Chiang, S. N. Yang, L. Tiator, M. Vanderhaeghen, 

and D. Drechsel, Phys. Rev. C68, 045202 (2003), nucl- 

th/0212106. 

[61] F. Tabakin, S. A. Dytman, and A. S. Rosenthal (1988), 
presented at Topical Conf. on Excited Baryons Troy, 
N.Y., Aug 4-6. 



The SAID interactive suite of codes is avail- [63] S. L. Adler, Phys. Rev. 137, B1022 (1965). 
able by secure shell interface via: ssh -C -X 
saidiSsaid . phys . gwu . edu. 



